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Abstract 

We revisit the classic problem of the stability of drops and jets held by surface tension, 
while regarding the compressibility of bulk fluids and spatial dimensions as free param- 
eters. By mode analysis, it is shown that there exists a critical compressibility above 
which the drops (and disks) become unstable for a spherical perturbation. For a given 
value of compressibility (and those of the surface tension and density at the equilibrium), 
this instability criterion provides a minimal radius below which the drop cannot be a 
stable equilibrium. According to the existence of the above unstable mode of drop, which 
corresponds to a homogeneous perturbation of cylindrical jet, the dispersion relation of 
Rayleigh-Plateau instability for cylinders drastically changes. In particular, we identify 
another critical compressibility above which the homogeneous unstable mode is predom- 
inant. The analysis is done for non-relativistic and relativistic perfect fluids, of which 
self- gravity is ignored. 
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1 Introduction 

The free oscillations of liquid droplets were studied by Kelvin [1] and Rayleigh [2] more than 
hundred years ago. Later, Lamb [3], Chandrasekhar [I], Reid [5] and others generalized the 
analysis to take into account the effects of an outer fluid, viscosity, and so on (see, e.g., Intro- 
duction of reference [6] for a brief but nice review and a more complete list of references). In 
the simplest case of an inviscid droplet in vacuum (or approximately in air), the droplet with 
unperturbed radius r and constant density p oscillates with the angular frequency given by 

Q 2 = ^(£- !)£{£+ 2), (1.1) 
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where t = 0, 1, 2, . . . denotes the mode number and a is the surface tension responsible for the 
oscillation. 

Another important phenomenon associated with the surface tension is the drop formation 
resulting from an instability of cylindrical jets. Theoretical studies date back to the early in- 
vestigations by Plateau [7], Rayleigh [8], and later by Chandrasekhar [9] (see [10] for a compre- 
hensive review). For example, an inviscid cylinder with non-perturbed radius r, the sinusoidal 
perturbation with wavenumber k evolves in time with the growth rate given by 



where I is the modified Bessel function of the first kind. This dispersion relation tells us that 
any cylinder that is longer than 2irr is unstable and the most unstable mode, which roughly 
determines the size of droplets forming, appears at wavelength A ~ 9r. 

The above two phenomena, the oscillations of droplets and the instability of cylinders, are 
not only of fundamental importance in the theoretical fluid mechanics but also important from 
industrial points of view. Therefore, they have been studied theoretically and experimentally 
in variety of physical situations considerable. However, the effect of non-zero compressibility or 
finite sound velocity of fluids has not been well studied, as far as the present author knows. The 
reason would be that the compressibility of liquids in many non-extremal situations is expected 
to be negligible and give rise no perceivable effects. 

Notwithstanding the above general expectation, in this paper we re-investigate the two 
classic problems, while allowing the bulk fluids to have non-zero compressibility or finite sound 
velocity. Our analysis reproduces the results (11.11) and (11.21) in the incompressible limit. The 
results are a little surprising and intriguing, at least, at theoretical level. The stability structure 
of droplets and cylinders with finite compressibility is rather richer than expected. It will be 
shown that there exists a critical compressibility above which the droplet (or a disk in a two- 
dimensional case) becomes unstable for a spherically symmetric perturbation. According to 
the existence of such an instability of disk, a cylinder, whose cross section is the disk, becomes 
unstable above the critical compressibility. These instability criteria can be interpreted as 
follow. For given parameters of the fluid and surface, namely the sound velocity c\ = dp/ dp, 
surface tension, and the density at the equilibrium, there exists a minimum radius of droplet 
and cylinder below which they cannot be stable equilibria!*]. Such a minimum radius is identified 

1 If one is familiar with theories of stellar structure, this result may not be totally surprising. That is, a 
self- gravitating fluid ball is unstable if the heat capacity ratio (or adiabatic index), denoted by 7 conventionally, 
exceeds 4/3 in Newtonian gravitational theory (see, e.g., QT]). This critical value is changed to a higher value 
in general relativity (see, e.g., [H]). Furthermore, if 7 is slightly larger than 4/3, then the spherical fluid ball is 



a kr (1 — (kr) 



2 )h(kr) 
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where c 2 s = dp/de in this case. 

Fortunately (or to the author's regret), in many non-extremal systems such as the water 
in air at room temperature r m j n is extremely small (r m j n ~ 10~ n m), and the instability 



however, the existence of the minimal radius is significant to prove that the Euler equation 
supplemented by Young-Laplace's stress balance relation at the surface, which governs the 
perfect-fluid systems considered in this paper, is not well defined for arbitrary values of (a, c s , p), 
while these parameters are usually supposed to take any positive finite values. 

The organization of this paper is as follow. In section [2] we review the Euler equation and 
Young-Laplace relation for non-relativistic fluids, and then we derive the dispersion relation for 
the perturbations of spherical and cylindrical equilibria. In section [3j we analyze the dispersion 
relation to obtain the stability criteria for the droplets (section 13.1 p and cylinders (section l3.2p . 
We generalize the analysis to the relativistic fluids in sections H] and [5j Section [6] is devoted to 
summary and discussion. The stability of the droplets for non-spherical (£ ^ 0) perturbations 
are shown in appendix [B] 

unstable to the radial perturbation provided the radius is less than a critical radius, which is much larger the 
Schwarzschild radius j!2| . The author thanks an anonymous referee for pointing out these points. 

2 We will consider axially symmetric fluids in arbitrary dimensions, by leaving n as a free parameter. The 
reason to do so is twofold. By leaving the spatial dimension as a free parameter, one can treat the disks, drops, 
and cylinders at once. Another reason is related to the so-called fluid/gravity correspondence [13) (see also [14j 
for a review and a complete list of references), that was found in the context of the string/M theories and 
relates the fluid mechanics in d dimensions to a gravitational theory in (d+ 1) dimensions. In this context, the 
spatial dimension often has to be treated as a free parameter and some phenomena such as the Raylcigh-Platcau 
instability and its subsequent dynamics have been known to depend crucially on the dimension [1 51 116) . It is 
added, however, that within the analysis in this paper, we could not find any qualitative difference originating 
from the difference of dimension. 

3 The author confesses that he cannot say for certain whether or not there are systems where r min is macro- 
scopic. See appendix |A"1 for discussions on the values of r m i n . 




From a theoretical or mathematical point of view, 
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Figure 1: An axially symmetric fluid and the cylindrical coordinates in a (n + 3) -dimensional 
spacetime M 1 '"-" 1 " 2 . The fluid surface at z = Const is a n-sphere S n . 



2 Non-relativistic perfect fluid with boundary 
2.1 Euler equation and Young-Laplace relation 

We consider compressible inviscid fluids in d-dimensional (d > 3) flat spacetime M 1 ^ -1 . Denot- 
ing density, pressure, and velocity field by p, p, and v 1 (I, J — 1, 2, . . . , d — 1), respectively, the 
continuity equation and Euler equations p2] are 

dtp + Vj(pu J ) = 0, (2.1) 
dtipv 1 ) + Vj(/w V + pg IJ ) = 0, (2.2) 



where V/ is the covariant derivative compatible with a flat metric gu |18j . 

We assume that a lump of fluid is supported by a constant surface tension a (> 0). Then, 
the Young-Laplace relation [17], describing the normal-stress balance at the surface, is given 
by 

P = aK\ f=Q . (2.3) 

Here, / is a scalar function with which the surface is identified by / = 0, and k is \{d — 2)-times] 
the mean curvature of the surface, given as the divergence of the unit normal vector n 1 (or as 
the trace of extrinsic curvature), 

Furthermore, we assume that the surface is convected with the fluid, which is expressed as 

(d t + v I V I )f = 0\ f=0 . (2.5) 

Since we are interested in axially symmetric fluids, it is convenient to work in cylindrical 
coordinates. Writing the spacetime dimensions as d = n + 3 (n = 1,2,...), the line element is 
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written a^| 



giJ dx l dx 3 = dz 2 + dr 2 + r 2 ds 2 n = S ab dx a dx b + r 2 ^{d)dB l dB\ (2.6) 



where x a := (z,r) and r y i j(6)d6' l d6 : > — 1, 2, . . . , n) is the line element on the unit n-sphere. 
Now, we assume that the fluid and its surface have SO(n+l) symmetry around the ^-coordinate. 
Then, the coordinate dependence of fluid quantities are given by 

p = p(t,z,r) v a = v a {t,z,r), v l = 0, f(t,z,r) = r-R{t,z), (2.7) 

where R(t, z) is the function representing the local radius of the fluid surface (see figured]). 
With the above ansatz, continuity equation (12. ip and Euler equation (12. 2p are written down 

as 

(d t + v z d z + v r d r )p + p (d z v z + d r v r + -v^j = 0, (2.8) 

p(d t + v z d z + v r d r )v a = -d a p. (2.9) 
The mean curvature of the surface, appearing in Young-Laplace relation (I2.3p . is given by 

K = n d " R ( 2 10) 

R[l + {d.R) 2 ] 1 / 2 [l + {d z R) 2 f/ 2 ' y } 

Finally, the kinematic boundary condition (I2.5P is reduced to 

d t R + v z d z R = v r \ r=R . (2.11) 

2.2 Linear perturbation of cylinders 

The above system obviously allows the static cylinder as an equilibrium solution, where the 
constant pressure po and the radius of cylinder R = r satisfy 

Po = a-. (2.12) 

Now, we perturb this equilibrium solution. We can assume that the perturbation results 
from a sinusoidal disturbance of the local radius given by 

R(t, z) = r [1 + ee" 1 cos(kz)] , (2.13) 

where |e| 1 is a small parameter. Such a disturbance leads to those of pressure and velocity, 
where the disturbed pressure at 0(e) may take form of 

p(t, z, r) = p [1 + ee^P^r) cos(kz)] . (2.14) 



4 For d = 3 we just discard the z-coordinate, and the cylindrical coordinates reduced to polar coordinates. 
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In general, from the perturbations of continuity equation (12. ip and Euler equation (I2.2p . the 
pressure perturbation has to satisfy the wave equation 



{d 2 - c 2 s V J Vj)8p = 0, (2.15) 

where c 2 := dpo/dpo (po is the density at the equilibrium and < c s < oo) is the sound velocity 
square of the bulk fluid, and we denote 0(£)-perturbation of any quantity X by 5X hereafter. 
Plugging expression f |2. 14j) into equation (I2.15p . we obtain 

* P + ^-(e + 4V = 0. (2,6) 



dr 2 r dr \ c 2 / 

With the regularity at the axis (r = 0), this equation is solved by the modified Bessel function 
of the first kind 

^Hcfc^>, + (2.17) 



* 



where C is an integration constant. The perturbation of Young-Laplace relation ( 12. 3p . 5p = 
adK\ r= ji, fixes the integration constant, 

C = J n -} kr ^fl~ l) '\ (2.18) 
nl(n-\)/2 {Kr ) 

The perturbation of Euler equation in the r-direction (12. 9p is pod t 5v r = —d r 5p. On the 
other hand, the velocity in the r-direction at the surface is given by 8v r = dt5R\ r= n from 
kinetic boundary condition ( 12. lip . The combination of these two yields 

d r 5p= - Po d 2 5R\ r=R . (2.19) 

Plugging equations ( I2.13p . ( I2.14p . (I2.17p . and (I2.18P into (12. 19j) . and eliminating a derivative of 
the modified Bessel function, we finally obtain the dispersion relation of perturbations for the 
compressible cylinder 

2 ° r n \2tt^ ^(n+l)/2 (-^ r o) fn 

u 2 = 3 [n - (kr ) 2 }Kr -± " ) . 2.20 

In the incompressible limit (c s — > oo), this dispersion relation reduces to equation (II. 2ft for 
n — 1 and the one derived in [19] for general n. 

To simplify the dispersion relation, we introduce here the following dimensionless quantities 

k := r k, u := ' u, (3:= f-^j ' c; 1 (> 0). (2.21) 
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This /3, being basically the reciprocal of the sound velocity, serves as the parameter representing 
the (adiabatic) compressibility of the bulk fluid. In terms of these dimensionless quantities, 
dispersion relation (I2.20P is equivalent to the relation between k and Cj for which the following 
function vanishes, 

F(Cj, k) : = Cj 2 - (n - P)K /(n+1)/2 W , K := (k 2 + feu 2 ) ^ . (2.22) 

3 Analysis of dispersion relation for the non-relativistic 
fluid 

We have derived dispersion relation (I2.20p for the perturbation of cylinders in M 1,n+2 . However, 
since the cross section of a cylinder is a sphere, equation (I2.20p with k — provides the 
oscillation frequency (or the growth rate of a possible instability) of droplets in M 1,n+1 (i.e., 
without the ^-direction) . Thus, we divide this section into two parts, subsections 13.11 and 
13.21 In the former, we investigate the dispersion relation with k — 0, which corresponds to 
both a spherically symmetric perturbation of droplets and homogeneous (in the ^-direction) 
perturbation of cylinders. Then, in the latter we investigate the dispersion relation for general 
modes with k > 0. 

3.1 Instability of drops (k = mode) 

Setting k = in fl2T20l) . we have 

F(u,0)=u -nfiu-z 7— rr. (3.1) 

J(n-1)/2(PW) 

It is noted that in the incompressible limit ((3 — > 0), the zero of F(u,0) identically vanishes 
(Cj = 0). This just says that the incompressible droplets cannot oscillate nor collapse (nor 
expand) with keeping the spherical symmetry. 

Let us see the behavior of F(ou, 0) for small Cj. Expanding equation f )3.ip around Cj = 0, one 

has 

n o2\ - 2 i n n4~ 4 



F(Cj, 0) = 1 - —B> Cj 2 + /3^ 4 + 0(u b ). (3.2) 

\ n + 1 / {n + iy{n + 3) 

From this, one can see F(Cj,0)\^ =0 = d&F(Cj, 0)|^ =0 = 0. On the other hand, it is easy to 
see that Um^ooF(Cj,0) = +oo from equation (13. ip . Therefore, if 8^F(Cj, 0)|& =0 < holds, 



S 



(a) (b) 




P r/r 



Figure 2: (a) Compressibility parameter (3 vs. dimensionless growth rate Cj of the radial pertur- 
bation of non-relativistic droplets in M 1,n+1 for several n. n — 1, 2, 3, 4 from the bottom to the 
top. (b) The radial eigenfunction P(r) in the n = 2 case for several values of compressibility 
parameter. /3//3 Cj i = 1.04, 1.13, 1.30, 3.00 from the bottom curve to the top. Negative P(r) 
corresponds to the decrease in the pressure, resulting from the increase in the radius of droplet. 



F(Cj, 0) must have at least one positive zero from continuity. From equation (13. 2p it is clear 
that <9|F(u>, 0)1^=0 < holds if (3 is larger than a critical value, 



n 1 



n 



& > Pa := • (3.3) 



In this case, from equation ( 13. 2 p the behavior of Cj near /3 c l is 

\2\ 1/4 



^±^1) (/3 _ p + o - /3 Cil ) 3/2 ) . (3.4) 
n + 1 / 



The global /3-dependence of cD obtained from equation (13. ip is shown in figure |2Ja)- One 
can see that in the large-/? limit, Cj increases linearly with /3. This asymptotic behavior, that 
is in fact Cj ~ n/3, can be derived from equation (13 .ip with using a property of modified Bessel 
function, lim^ 00 I( n+1 y2(PCj)/I( n _ 1 y 2 ((3Cj) = 1. Radial function P(r) in the n = 2 case for 
several values of /3 is shown in figure Mjo). One can observe the non- uniform (in the r-direction) 
decrease in the pressure, resulting from the increase in the droplet radius. The non-uniformity 
of the pressure perturbation is amplified as the compressibility increases, which is interpreted 
that the compressibility reduces the propagation speed of density fluctuations generated near 
the surface. 

With using equations (I2.12p and 02.2ip . the instability criterion (I3.3P can be written in 
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ro < r min := ; (3.5) 



several forms. Namely, c 2 s < po/[(n + l)po] or equivalently 

n a 
a + 1 p cj 

Inequality (13. 5 p is striking, which means that there exists a minimum radius r min only above 
which the droplet and cylinder can exist stably. The critical radius is determined by three 
parameters (po,a,c s ), which do not restrict each other at least from a macroscopic point of 
view, although should be correlated microscopicalljjE 

One could consider non-spherical perturbations of the drop by allowing ^-dependence of R 
in equation (12.131) . As in the incompressible case, however, the non-spherical modes turn out 
to be oscillatory for the compressible fluids too, proving that the droplets are stable for such 
perturbations. See appendix IB. II for a proof. 



3.2 Rayleigh-Plateau instability (k > modes) 

We proceed the dispersion relation (12.201) for general values of k > 0. Since the homogeneous 
(k = 0) mode becomes unstable (u > 0) above the critical compressibility f3 Cj i, one can expect 
the behavior of the dispersion relation for k > changes at the critical compressibility. 

The dispersion relation (12.201) for n = 1 is shown in figure [3] for several values of /3. Quali- 
tative behaviors are independent of n. For < f3 < /3 C) i, the dispersion relation is qualitatively 
the same as the usual dispersion relation of Rayleigh-Plateau instability (for an incompress- 
ible inviscid fluid), although the growth rate for all k G (0, y/n) mildly increases with 0. The 
wavenumber and growth rate of the most unstable mode at = /3 c l obtained numerically are 
given in table [TJ As (3 exceeds /? C) i, the homogeneous (k = 0) mode begins to have positive 
growth rate as shown in section 13.11 Incidentally, the dispersion relation deviates from that 
of usual Rayleigh-Plateau instability. Note that even if /3 exceeds /3 Cj i only slightly, the most 
unstable mode is still in k > 0. However, there exists another critical value of ft [we call it 
(3 C ,2 (> 0c,i)] above which the homogeneous mode becomes the most unstable one. The values 
of /3 Ci 2 and growth rate of the most unstable mode therein, that we denote by w* |2) are also 
given in table [TJ They can be obtained by solving numerically the following coupled algebraic 
equations for (3 and &, 

F(u;,k)\ k=o = 0, dfF(£,% =0 = 0. (3.6) 

5 The value of r m in for the water in air at temperature 25 °C is around 2.13 x 10~ n m (estimated with 
n = 2, a = 72.0 x 10~ 3 J/m 2 . p = 1.00 x 10 3 kg/m 3 , and c s = 1.50 x 10 3 m/s [20]), where the fluid-mechanical 
description has already broken down. Thus, the instability found plays no central role for such a fluid. It would 
be interesting to look for a system in which r m i n is larger than or comparable with the length scale where the 
fluid approximation breaks down. See also appendix [X] 
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We note that /3 C .2 is of course larger than /? C; i but only slightly. Thus, we can say that for 
generic values of (3 larger than /3 c l the homogeneous unstable mode is predominant. 

0.8 
0.6 
<3 0.4 
0.2 

°0 0.2 0.4 0.6 0.8 1 

k 

Figure 3: Dimensionless wavenumber k vs. dimensionless growth rate Cj of the instability of 
non-relativistic cylinder in K 1 ' 3 (n = 1) for several values of compressibility. = 0.00, 1.41 (= 
/3 Cj i), 1.45, 1.49 (= /3 C)2 ), 1.53 from the bottom to the top. 




Table 1: The critical values of compressibility parameter, /3 c l and /3 C) 2, and characteristic 
dimensionless wavenumber and growth rate. /3 Cj i is the compressibility above which the cylinder 
is unstable for the homogeneous (k = 0) perturbation. and are the wavenumber and 
growth rate of the most unstable mode when f3 = /3 Cj i, respectively. (3 Ct2 is the compressibility 
above which the homogeneous (k = 0) mode becomes the most unstable one. u)* ; 2 is the largest 
growth rate at (3 = /3 C)2 . 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




1.414 


1.225 


1.155 


1.118 


1.095 


1.080 


1.069 


1.061 


1.054 


1.049 


k*,i 


.4388 


.6252 


.7596 


.8672 


.9580 


1.037 


1.107 


1.171 


1.229 


1.282 




.5244 


.8504 


1.115 


1.345 


1.552 


1.742 


1.919 


2.086 


2.244 


2.394 


Pc,2 


1.489 


1.293 


1.215 


1.172 


1.144 


1.124 


1.110 


1.098 


1.089 


1.081 


&*,2 


.6388 


1.039 


1.358 


1.632 


1.875 


2.097 


2.301 


2.492 


2.671 


2.840 
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4 Relativistic perfect fluid with boundary 



The argument in the preceding sections can be generalized to relativistic fluids. We consider 
relativistic fluids in the d- dimensional flat spacetime IR 1 '^ -1 (d > 3) with the spacetime coordi- 
nates = (ct, x 1 ) (/i, v — 0, 1, 2, . . . , d — 1; I, J — 1, 2, . . . , d — 1), and denoting the flat metric 
by [with the so-called almost-plus notation (— , +, • • • , +)]. The speed of light is set to unity 
(c = 1). Symbols appearing hereafter have the same meanings as in the non-relativistic case, 
otherwise noted. 

4.1 Euler equation and Young-Laplace relation 

The energy-momentum tensor of a relativistic perfect fluid held by surface tension (see, e.g., 
[211122]) is given by 

= (ewV + pP^)e(-f) - a\Vf\h" v 5{}). (4.1) 

Here, e is the energy density; the normalized G?-velocity field (u^u^ = —1); the Heaviside 
step function; S the Delta function; P^ v := g^ u + u^u u and W := g^ v — n^n u are the projection 
tensors; is the unit normal of the surface defined by = V At //(V' y /V i/ /) 1 ' 2 ; and are 
orthogonal each other at the surface u^n^ = 0|/ =0 . 

Projecting the energy-momentum conservation V = onto and P^, we obtain the 
relativistic continuity and Euler equations, 

u a V a e + (e + p)V a u a = 0, (4.2) 
(e + p)u a V a u» = -P» a V a p. (4.3) 

The Young-Laplace relation, obtained by projecting the surface contribution of the energy- 
momentum conservation onto n M , takes the same form as the non-relativistic case ( 12. 3p . but 
the mean curvature in the this case is given by the d- dimensional divergence, 

k = V/. (4.4) 

The kinematic boundary condition is given by 

u a V a f = | /=Q . (4.5) 

As in the non-relativistic case, we introduce cylindrical coordinates in which the line element 
of flat spacetime is given b}|^| 

g^dafda? = -dt 2 + dz 2 + dr 2 + r 2 ds 2 n = rj ab dx a dx b + r 2 j ij (9)d9 i d9 j , (4.6) 
6 For d = 3we just discard the z-coordinate as in the non-relativistic argument. 
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where x a = (t,z,r) and r] a b = Diag. (— 1, 1, 1) is the three-dimensional flat Lorentzian metric. 
Assuming that the fluid and surface are axially symmetric around the z-axis, the coordinate 
dependence of the pressure, velocity field, and surface are 

p = p{t,z,r), u a = u a {t,z,r), u* = 0, f = r-R(t,z). (4.7) 

Continuity equation ( 14. 2ft and Euler equation (14. 3 p are written down aJll 

u a d a e + (e + p) (d a u a + V) = 0, (4.8) 
(e + p)u b d b u a = -P*d b p, (4.9) 

where the indices (a, b, . . .) are raised and lowered by ?7 afe and 77^, respectively. The mean 
curvature ( I4.4p . appearing in the Young-Laplace relation (I2.3p . is given by 

n [1 - (d t R) 2 )d 2 z R - [1 + (d z R) 2 }d 2 R + 2(d t R)(d t d z R)d z R 



R[l - {d t R) 2 + {d z R) 2 Yl 2 [1 - (cW + {d z R) 2 f/ 2 

(4.10) 

Finally, kinetic boundary condition (14. 5p reads 

u%R + u z d z R = u r \ r=R . (4.11) 

4.2 Linear perturbation of cylinders 

In general, perturbing relativistic continuity equation (14. 2p and Euler equation H4.3[) around a 
static equilibrium where (p, e, u^) = (p , e , 6^), we have 

d t 5p + c 2 s (e +p )VJu a = 0, (4.12) 
(e„ + Po )d t 5u» + PffiVJp = 0, (4.13) 

where P^ := g^ + and c 2 := dp /de is the sound velocity square (note that c s < 1 from 
the causality). Eliminating 5m m from these two equations, one obtains an wave equation for the 
pressure perturbation, 

- <?.P$V a Vp) 5p = 0. (4.14) 

As in the non-relativistic case in section 12.21 the axially symmetric relativistic system in 
section 14.11 allows the cylinder as a static equilibrium, where constant pressure po and radius 



' In the present coordinates the non- vanishing components of Christoffel symbol are r°- = —rS^jij, Vj a = 
r~ 1 Sj5a, and TJ fc = ^ 7 ^rj fc , where ™rj fc is the Christoffel symbol with respect to 7^. With using these, the 
vector derivatives are calculated. For example, V a u& = <9 a it{,, V^tij = ru r jij, and X7 a u a = d a u a + nr -1 ^ . 
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R = r satisfy equation (I2.12p . Plugging ansatz f |2.14|) into wave equation (I4.14p and using 
the perturbation of Young-Laplace relation 5p = a5K\ r= R, we obtain the radial function in the 
relativistic case, 

f(r)= >-^y j ^ r » , K:=(e + $y. (4,5) 

The perturbation of Euler equation (I4.13P in the r-direction reads (e +po)d t Su r + d r 5p = 0. 
On the other hand, the perturbation of kinetic boundary condition (14. lip reads dtSR = 5u r \ r= R. 
Eliminating 5u r from these two equations, we obtain 

(eo+p )%6R=-d T 6p\ r=R . (4.16) 

Plugging equations (12.131) . (12.14[) . and (I4.15|) . into (I4.16p . we obtain the dispersion relation for 
the perturbation of relativistic compressible cylinder 



u = —— — —^[n-ik +u)r \Kr - — — . 4.17 



It is noted that a similar result for a relativistic compressible fluid with a particular equation 
of state was obtained in |16j . 

We introduce the following dimensionless quantities (remember that we have already set 
the speed of light to unity) , 

k;=r k, u:=r u, a := (ebr )"V, (3 := c; 1 (> 1). (4.18) 

Then, dispersion relation (14.171) is the relation between k and Cj for which the following function 
vanishes, 



F(u,k) :=uj 2 ^[n-(k 2 + u 2 ))K I(n+1)/2iK \ K := (k 2 + f3 2 Cj 2 ] 1 ' 2 . (4.19) 



J {n+1)/2 (K) fs._(l2 , Q 2,~.2\ l/2 
l(n-l)/2{K) 

Function F(u, k) depends on not only the dimensionless compressibility (3 but also the di- 
mensionless surface tension a. This is contrast to the non- relativistic counterpart (I2.22p . that 
depends only on the compressibility (3. 



5 Analysis of dispersion relation for the relativistic fluid 

As explained at the beginning of section I3.1[ the homogeneous (k = 0) mode has the special 
meaning that it corresponds to both the homogeneous perturbation of the cylinders in M 1,n+2 
and the spherical perturbation of droplets in R 1,n+1 . In the first subsection below we look into 
the k = mode, then we proceed the analysis of general k > perturbations in the second 
subsection. 
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5.1 Instability of drops (k = mode) 

Setting k = in equation (14.191) . we have 



-./(n+l)/2 



F(o>, 0) = & + — — (tf - n)pCj '™'J-' (5.1) 

Let us see the behavior of F(u,0) for small Cj by expanding it around Cj — 0, 

^ 2 \ a2 n/3 2 + (n + l)(n + 3) „„ 2 „ 4 



0) = 1 - t 7~\~7i tt/3 2 Cj 1 + 7 ' , ' ' a/ny + 0(w b ). (5.2) 

v 7 V (n+l)(l + no-) y (n+ l) 2 (n + 3)(l +na) \ j \ j 

From this, one can see F(Cj,0)\& =0 = d^F(u,0)\^ =0 = 0. On the other hand, one can see 
lirriu^coF {fa , 0) = +oo from equation (I5.ll) . Thus, if d?F(u, 0) |^=o < holds, F(u>, 0) must have 
at least one positive zero from continuity. From equation (15.2ft . one can see that &$F(Cj, 0) \&=q < 
holds if the compressibility parameter is greater than a critical value, 

P>^:=( {n+1)(1 . +na) Y 2 . (5.3) 



na 

In this case, the /3-dependence of growth rate near (3 = (3 Ct i can be read from equation (15. 2p . 

» = ( (n + l)(l+'yr£, + 3K]j ^ - A '' ),/2 + ° «" - ^) • (5-4) 
The global behavior of Cj((3) for several values of a is shown in figure H^a). In the present rela- 
tivists case, in contrast to the non-relativistic case, the growth rate asymptotes to a constant 
in the large-/? limit, that is Cj ~ y/n in fact. This can be derived from equation ( 14. 191) . Radial 
function P(r) is shown in figure HJ^b). 

Getting back to dimensionful quantities with equation (I4.18p . instability criterion ( 15. 3 p is 
rewritten as 

r < r min := n a (1 - (n + l)c 2 s ) . (5.5) 
n + 1 e c| 

Namely, there exists minimum radius r min below which the drops and cylinders become unsta- 
ble. See appendix IB. 21 for the proof of stability for non-spherical perturbations. 

5.2 Rayleigh-Plateau instability (k > modes) 

The behavior of dispersion relation cD(fc) is quite similar to the non-relativistic case, except that 
the critical values of compressibility, /3 c l and /3 C)2 , depend on a. Numerical plot of Cj — d(k) 
for several values of (3 is shown in figure EJa). The value of a does not affect the qualitative 
behavior of the dispersion relation. In order to see that the second critical compressibility /3 C> 2, 
above which the k = mode is the most unstable one, is only slightly larger than the first 
critical value (3 C) i in all range of a, we numerically plot j3 C) i and (3 Ct2 for n — 1 in figure [3(b). 
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(b) 




fi r/r 

Figure 4: (a) Compressibility parameter /3 vs. dimensionless growth rate of the radial pertur- 
bation of droplet in R 1,3 (n = 2) for several values of dimensionless tension, a = 0.100 (thick 
solid) 0.300 (thick dashed), 1.00 (thin solid), 20.0 (thin dashed), (b) Radial function P(r) for 
the same set of a as that in figure H](a). n = 2 and /3//3 c ,i = 1-50. 



(b) 




Figure 5: (a) Dimensionless wavenumber k vs. dimensionless growth rate u of the perturbation 
of relativistic cylinder in M 1,3 (n = 1) for several values of compressibility. (5 = 1.00, 2.00 (= 
Pc,i), 2.13, 2.27 (= Pc,2), 2.45 from the bottom to the top. The dimensionless surface tension 
is fixed to a = 1. (b) The a-dependence of the critical compressibility, /3 Ci i (solid) and (3 Ct 2 
(dashed), for n — 1. 
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6 Conclusion 



We have investigated the stability of spherical drops and cylindrical jets held by the surface 
tension, in particular, its dependence on the compressibility f3 or sound velocity c s of the 
bulk fluids. For simplicity, we have focused on the perfect fluids (i.e., inviscid fluids with no 
heat transfer) immersed in vacuum, while we consider both the non-relativistic and relativistic 
fluids in general dimensions, which allow us to treat the disk, droplets, and cylindrical jets 
simultaneously in a systematic way 

As the main result, we have shown that there exists the critical compressibility /3 Cj x for 
both the non-relativistic and relativistic fluids, equations (13. 3 j) and (15. 3p . above which the 
spherical drops are unstable for the spherical perturbation. For given parameters of the fluid 
and surface, i.e., the surface tension a, and the sound velocity c s and density at the equilibrium 
Pq (or e in the relativistic case), the instability criteria poses the lower limit on the droplet 
size r m i n ~ & / (.Po c l) [ see equations (13.51) and (15.51) ]. below which any droplets cannot be stable 
equilibria. 

We have shown also that according to the instability of disks and droplets, which corre- 
sponds to the instability of cylinders for homogeneous perturbations, the dispersion relation of 
Rayleigh-Plateau instability exhibits the significant change. Namely, for (3 > /3 c l the cylinders 
are unstable for the perturbations that are homogeneous in the axial direction. Furthermore, 
such a mode becomes the most unstable one above the second critical compressibility /3 C) 2, which 
is slightly larger than /3 c> i in general. 

Here, let us stress the significance of minimum radius r min . In the framework of fluid 
mechanics, it has been assumed or simply believed that any positive finite values can be given 
to the three quantities a, c s , and po (although they should be correlated each other if one 
pursues their origins from a microscopic point of view). We have shown, however, that the 
spherical droplets, which are the most fundamental equilibrium states of localized fluids in 
fact, exhibit the instability for r < r min ~ cr/(p c^). Therefore, one cannot give values to the 
three parameters freely in order to describe arbitrarily small-scale dynamics successfully. In 
other words, the systems defined by the Euler equation and Young-Laplace relation intrinsically 
contain the instability, and are not well defined in certain regimes of parameter space. 

We have adopted several assumptions for simplicity, such as the absence of viscosities, 
heat transfer, and outer fluids, the constancy of surface tension, and so on. In addition, the 
instability discovered is just the result of mode analysis, that can never predict the following 
dynamics. Therefore, there are many directions to proceed by generalizing the analysis in this 
paper. It would be interesting to see how the viscosities affect the instability. The nonlinear 
dynamics would be interesting too, even within the perfect-fluid approximation. 
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A Comments on r 



mm 



After submitting a draft of this paper, anonymous referees gave the author several useful 
comments and suggestions about the values of r min in actual physical systems. Here, some of 
these suggestions are noted for further studies. 



General order estimate. On dimensional basis, one can argue that the critical radius r m j n 
is of the order of a microscopic length scale as follow. Let us consider the relativistic case for 
simplicity. In many physical systems, it would be possible to assume c 2 s = dpo/de® ~ 1, which 
is equivalent to assume that the compressibility is proportional to the inverse of the energy 
density, e^deo/dpo ~ . In this case, r min in equation (15. 5p is reduced to r min ~ cr/e . Here, 
let us assume further that the surface tension is microscopically proportional to the energy 
density in the bulk as a ~ Ze , where / is a microscopic length scale. Such a length scale I could 
be a mean inter- molecular distance^. With these assumptions, one obtains r m j„ ~ /, which 
suggests that the instability is irrelevant since the hydrodynamic description itself breaks down 
at the scale of I. 



Liquid-drop model of nucleus. It is widely known that many features of nucleus such as the 
global behavior of binding energy, surface oscillations, and nuclear fission, can be understood 
with the liquid-drop model, in which the nucleus is modeled by a liquid drop of an incompressible 
(at leading order) fluid [24"l 125]. However, the compressibility, which allows radial oscillations 
of a drop {i.e., the so-called breathing mode), is important since it is directly related to the 
equation of state of nuclear matter, and necessary to the accurate estimate of nuclear properties 
(radii, masses, giant resonances, etc.). Though the argument in the above paragraph seems to 
show the instability found in this paper is irrelevant to nuclei, it would be interesting to compare 
systematically the parameters in the liquid-drop model and those in this paper. Incidentally, 
the confine/deconfine phase transition in the quantum chromodynamics (QCD) is expected to 

8 This should be justified in a molecular theory of capillary force. See, e.g., [25] . 
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be of the first order, and so the deconfined phase, i.e., the quark-gluon plasma (QGP), could 
exist as a drop of fluid around the critical temperature. Thus, it would be also interesting to 
consider the effect of compressibility to the QGP balls. 

Granular matter. Recently, it was reported that a kind of granular matters such as glass 
beads of tiny radius exhibit effective compressibility [26] and surface tension |27^l . The dynam- 
ics of granular matters cannot be described by hydrodynamics in general, and furthermore the 
origin of the surface tension is different from that in fluids. Therefore, one cannot apply the 
result in this paper as is to granular matters. It would be interesting, however, to examine the 
possibility that the new capillary instability found in this paper or its analog is effective or not 
in granular matters. As far as the author knows, the sound velocity of granular matter could 
be relatively small (e.g., c s ~ 1 m/s in an experiment [28]). Thus, there remains the possibility 
that the r min oc c~ 2 could be large. 



B Stability of drops for non-spherical perturbations 
B.l Non-relativistic case 

We show the stability of non-relativistic droplets in M 1,n+1 (n > 1) for non-spherical per- 
turbations [and the stability of cylinder in M 1,n+2 for homogeneous (in the r-direction) but 
non-spherical perturbations] . 

We work in the polar coordinates in which the line element of the flat space is given by 

g IJ dx I dx J = dr 2 + r 2 lij (9)d9 l de j . (B.l) 

In this coordinates, continuity equation ( 12. ip and Euler equation ( 12. 2 p in the r- and #*-directions 
are 

0, (B.2) 

-d rP , (B.3) 

~^d jP , (B.4) 

9 The existence of surface tension (capillarity) in a granular matter might be surprising since the attractive 
force between grains is much smaller than other forces at play (gravity, friction, inelasticity). In fact, the surface 
tension of the granular matter in experiment [27] (grass beads) stems from not the attractive force between the 
beads but a strong interaction between the beads and surrounding air. 



(d t + v r d r + v%)p + p (d r v r + V + DiV* 

\ r 

p((d t + v r d r + v l di)v r — r'y ij v l v j 



p[d t + V r d r + - + vWj ) V 1 
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where Di is the covariant derivative compatible with 7^. If we parametrize the scalar function 
/ as 

f(t,r t e) = r-R(t,e 1 ,e 2 ,...,e n ), (b.5) 

the mean curvature, appearing in Young-Laplace relation ( 12. 3p . reads 

n R 2 D 2 R+(D 2 R-R)(DR) 2 -(D i R)(DiR)D i D j R 

K ~ R[l + R~ 2 (DR) 2 y/ 2 R*[l + R~ 2 {DR) 2 } 3 / 2 ' ^ ' 

Here, (DRf := ^{D^DjR and D 2 R := fW^R. Kinetic boundary condition is 

d t R + v i d i R=v r \ r=R . (B.7) 

Obviously, equations of motion (IB.2j) . (lB.3p . and (1B.4|) . and boundary conditions (12. 3p [with 
equation ( IB. 61) ] and ( IB. 71) allow the spherical drop as a static equilibrium, where the constant 
pressure p and radius of sphere R = r satisfy equation (I2.12p . Now, we consider linear 
perturbations of this equilibrium resulting from the disturbance of surface, 

R(t,0) = r Q [l + ee wt Y{6)}. (B.8) 

Here, Y(8i,8 2 , . . . , 9 n ) is the harmonic function on the unit n-sphere, 

[D 2 + £(£ + n- 1)]Y(0) = 0, £ = 0, 1, 2, ... . (B.9) 

The perturbed pressure to 0(e) may be written as 

p(t, r, 9) = po[l + se wt P(r)Y(9)\. (B.10) 

Substituting expression (IB. 101) into wave equation (12 . 1 5j) . we obtain 

gp + 2 g_^ + <(< +r iA P = ft (B . U) 

ar z r ar \Cg r z J 

With the regularity at the origin and the perturbed Young-Laplace formula 5p = a5n\ r= R, one 
finds the following solves equation (IB. lip . 

p(r) = ^-m + n)r { r l)l2 h +{ r^)/2^r/c s ) _ 
nh +{ n-i)/2 (ur /c s ) r^" 1 )/ 2 

Substituting equations (l2~T2|) . ( TB~8]) . (|B~T0|) . and ( TB~T2|) into ( 12~T9|) . we obtain 

2_ °~ 1V « , \ { o wr /< + ( n+ i)/ 2 (wro/c s 



a; 2 = -—J£ -1)(£ + n) I I + — — / \ ■ (B.13) 
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If one takes n = 2, c s — y oo, and writes to — y ifi (i = V - 1) i n this relation, one reproduces the 
classic formula (II. ip . i.e., the angular frequency of oscillations for the incompressible inviscid 
droplet immersed in the three-dimensional vacuum [2J. 

If equation ( 1B.13[) has a positive root for a given £ (> 0), the spherical droplet is unstable for 



the perturbation labeled by £. Introducing the dimensionless quantities as in equation (12.21 j) . 
the problem is equivalent to find positive zero of the following function, 

Ftp) := u 2 + (i- !)£{£ + n ) + (£-!)(£ + n) ^ n+1)/2 ^% (B.14) 

^+(n-l)/2(P^J 



Note that Fq(u) is nothing but F(6j,0) in section [37X1 which was shown to have the positive 
zero. Taking into account the positivity of the modified Bessel function, one can easily see 

F e (u) > for £>l and to > 0. (B.15) 

Thus, there is no positive zero of Fi(£j) for i > 1, proving the stability of droplet for the 
non-spherical perturbations. 

B.2 Relativist ic case 

We show the stability of relativistic droplets in R 1 '™ 4 " 1 (n > 1) for the non-spherical perturba- 
tions [and that of relativistic cylinders in M 1,ri+2 for non-spherical but homogeneous (k = 0) 
perturbations]. 

We write the flat metric as 

g^dx^dx u = -dt 2 + dr 2 + r 2 y lj d9 l d9 j . (B. 16) 

Then, we perturb the spherically symmetric static equilibrium, where the constant pressure 
p and radius R = r satisfy equation (12.121) . with the same ansatz (1B.8[) and (IB.lOp as the 
non-relativistic case. The mean curvature (14.40 reads 



+ D4ri , aD „, ■ D - 2 ^ 21 ,/ 2 \R 2 + (DR) 2 ]R 2 d 2 R 



R[l - (d t R) 2 + R~ 2 (DR) 2 } 1 / 2 R 4 [l - (d t R) 2 + R~ 2 (DR) 2 ] 3 / 2 

- [1 - (d t R) 2 }R 2 D 2 R + (R— D 2 R)(DR) 2 + (D i R)(D ] R)D i D j R) . (B.17) 



Substituting expression ( IB.lOp into wave equation (I4.14p . one finds that radial function P(r) 
satisfies the same form of equation as ( IB. lip . With the perturbed kinetic boundary condition, 
one can fix an integration constant to obtain 

p(r) = IK) 2 + (£-!)(£ + n)]rt m FHn-m(ur/c s ) 
nli +{n -i )/2 (ujr /c s ) r(™-!)/ 2 
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Substituting equations flB.8j) . (IB. 10|) . and flB.18j) into (]4.16j) . one obtains 



a 



ur h+(n+i)/2(ur /c s ) 



u, 2 = - r — — - 3 [(ur y + (£-l)(£ + n)][£ + Z^ 'J^M^Zl^Zl ) . (B .l9) 

In terms of the dimensionless quantities in equation (14.181) . to find a positive root of the 
above equation is equivalent to find a positive zero of the following function, 

F&) := co 2 + -f^ 2 + (£-!)(£ + n)) U + ^ f ( " +1)/2 |f| ) . (B.20) 



Observe that Fq{CS) is nothing but F(uj, 0) in section I5TTI which was shown to have the positive 
zero. One can easily show with using the positivity of the modified Bessel functions that 

Fi(u) > for £ > 1, a > 0, and u > 0. (B.21) 

Thus, Ft(u) for £ > 1 has no positive root, proving the stability of relativistic droplets for the 
non-spherical perturbations. 
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